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Abstract. We present an X-ray study of the galaxy group RGH 80, observed by XMM-Newton. The X-ray 
emission of the gas is detected out to ~ A62h^Q kpc, corresponding to ~ 0.45r2oo- The group is relatively gas rich 
and luminous with respect to its temperature of 1.01 ± 0.01 keV. Using the deprojected spectral analysis, we find 
that the temperature peaks at ~ 1.3 keV around 0.11r2oo, and then decreases inwards to 0.83 keV at the center 
and outwards to ~ 70% of the peak value at large radii. Within the central ~ 60 kpc of the group where the gas 
cooling time is less than the Hubble time, two-temperature model with temperatures of 0.82 and 1.51 keV and 
the Galactic absorption gives the best fit of the spectra, with ~ 20% volume occupied by the cool component. 
We also derive the gas entropy distribution, which is consistent with the prediction of cooling and/or internal 
heating models. Furthermore, the abundances of O, Mg, Si, S, and Fe decrease monotonically with radius. With 
the observed abundance ratio pattern, we estimate that ~ 85% or ~ 72% of the iron mass is contributed by SN 
la, depending on the adopted SN II models. 
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1. Introduction 

In the past decade X-ray imaging spectroscopy has be- 
come a powerful probe of the physical conditions of the 
hot (r^lO^ K) intergalactic medium (IGM) of galaxy 
groups. Spatial distributions of gas temperature, density 
and metallicity, as well as the structures of the gravita- 
tional potential and dark matter halos have been mea- 
sured for numerous galaxy groups. This improves signif- 
icantly our knowledge of dynamical properties and for- 
mation of galaxy groups. Among these, the discovery of 
similarity breaking of some IGM properties from clusters 
to groups is crucial for the study of cosmic evolution of the 
hot IGM in very massive halos of the universe. The well- 
known example is the deviation of the X-ray luminosity 
(Lx) and temperature (T) relation of groups and clus- 
ters from the prediction of self-similar model (e.g. Edge 
& Stewart 1991; David ct al. 1993; Wu, Xue & Fang 
1999; Helsdon & Ponman 2000; Xue & Wu 2000 and ref- 
erences therein). The break lends strong support to the 
argument that the IGM of groups was much more af- 
fected by non-gravitational processes, such as radiative 
cooling and/or feedback of star formation, than the IGM 
of clusters. Another example is the detection of entropy ex- 
cess in clusters and groups (Ponmam, Cannon & Navarro 
1999; Lloyd-Davies, Ponman & cannon 2000; Xu, Jin & 
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Wu 2001), which reinforces the presence of strong non- 
gravitational effects especially in the inner regions of poor 
clusters and groups. The reliable Lx — T relations and 
the entropy profiles obtained from spatially resolved X- 
ray imaging spectroscopy may allow one for the first time 
to distinguish between the competing non-gravitational 
models (Voit et al. 2002, 2003; Mushotzky et al. 2003). 

X-ray spectra of groups are dominated by emission line 
features. The wealth of characteristic emission lines pro- 
vides a robust tool for understanding the metal enrich- 
ment processes in the IGM and possibly for constraining 
the supernova history of galaxy groups (e.g. Finoguenov 
& Ponman 1999). However, an accurate determination of 
the metal abundances depends strongly on precise mea- 
surements of the temperature structure of the hot IGM in 
groups, as the dominated Fe L-shell lines are highly tem- 
perature sensitive (e.g. Buote 2002; Buote et al. 2003a). 
The high spectral resolution of XMM-Newton allows us 
to obtain the abundances of O and Mg in addition to Si, 
S and Fe, and provides us a deep insight into the chemi- 
cal evolution history of groups with the abundance ratio 
pattern and the metal mass to light ratios (e.g. Xu et al. 
2002; Matsushita, Finoguenov & Bohringer 2003). 

In this paper, we report the results from new XMM- 
Newton observations of RGH 80, which was first identified 
by Ramella, Geller & Huchra (1989) in the center for as- 
trophysics redshift survey. It was identified as an extended 
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X-ray source in the ROSAT All-Sky Survey with an X- 
ray luminosity of 3.33 x lO^^h^^ erg s'^ in the 0.1-2.4 
keV band within an aperture of l-Qh^^ Mpc (NRGs241; 
Mahdavi et al. 2000). Global fitting of the ASCA SIS 
and GIS data by Davis, Mulchaey & Mushotzky (1999) 
yielded a mean gas temperature of 1.02 ± 0.05 keV, an 
average abundance of 0.26^0 ^8^0 for the a elements and 
O.2Ol[J;[56^0 for iron. However, Buote (2000) pointed out 
that a two-temperature spectral model provides a better 
fit to the ASCA spectra of RGH 80 and has a metallic- 
ity that is substantially higher than that obtained by the 
single temperature spectral model, Z = 0.67^0^24-^0 • We 
will use the XMM-Newton data to examine the spatially 
resolved X-ray properties of the group and explore their 
physical applications. 

In section 2, we describe the observations and data re- 
duction procedure. We present the detailed spectral anal- 
ysis in section 3. In Section 4 and 5, we derive the gas 
density profiles, and calculate the three dimensional dis- 
tributions of gas and dark matter, respectively. We discuss 
and summarize the implications of our results in section 
6 and 7. Throughout this paper, we take = 50 km 
s~^ Mpc~^, r^M = 1 and f^A = 0. At the group redshift of 
0.037, 1' corresponds to 60.2 kpc. Unless stated otherwise, 
we adopt 68% confidence limit for error analysis. 

2. Observation and data preparation 

RGH 80 was observed with XMM-Newton for 33 ks in 
revolution 563 on January 5th, 2003. The two MOS-CCD 
cameras were operated in full frame mode, and the pn- 
CCD camera was operated in extended full frame mode. 
All cameras were covered with the Thinl filter. We gen- 
erated the calibrated events lists for the data by using 
the tasks emchain and epchain that are packaged in the 
XMM SAS v5.4.1 software. In the analysis, we keep the 
events with PATTERNS 0-12 for the MOS cameras, and 
the events with PATTERNS 4 for the pn camera. 

The EPIC background is mainly composed of solar 
soft protons, cosmic rays, and cosmic X-ray photon back- 
ground (Lumb et al. 2002; Read & Ponman T.J. 2003). 
The background induced by soft protons is time variable, 
and hence causes large variations of intensity (flares) in the 
light curves. It can be subtracted by discarding those flare 
periods. Since in the very high energy band the effective 
area of XMM is negligible and the emission is dominated 
by the particle background, we extract the light curves in 
the energy band of 10-12 keV and 10-15 keV for the MOS 
and pn data, respectively, in 100s bins. An inspection of 
the light curves docs not reveal any strong temporal vari- 
ations. In order to determine the mean count rate and its 
error (ct) for the quiescent periods, we recursively clean 
the light curve by rejecting those time intervals during 
which the counts are 3(7 outside the mean value, until the 
mean counts remain constant. We then define the thresh- 
olds at the ±3(7 level, and reject any time intervals outside 
these thresholds. After applying this screening criterion. 



the final useful exposures are 33.0 ks for MOSl, 32.7 ks 
for M0S2, and 26.6 ks for pn, respectively. 

We utilize the EPIC blank sky event files available 
from the XMM calibrations (Lumb 2002) to represent the 
remaining quiescent background, which is mainly com- 
posed of the cosmic X-ray background and the background 
induced by cosmic rays. We applied the same PATTERN 
selection and light curve screening criteria as have been 
used above to the blank sky events. Since some internal 
background fluorescence lines show strong spatial inhomo- 
geneities, such as Al K and Si K fluorescent lines in the 
MOS cameras (Lumb et al. 2002; Read & Ponman 2003) 
and Cu K fluorescent lines in the pn camera (Freyberg, 
Pfeffermann & Briel 2001; Lumb et al. 2002; Read & 
Ponman 2003), we use the sky cast script ^ to cast the 
background files into the corresponding sky coordinates of 
the source observations to ensure that the products of the 
background and the source are extracted from the same 
location. 

Before subtracting the background, we first correct the 
vignetting effect for both source and background data sets 
with the SAS task evigweight, which computes a weight 
coefficient for each photon by using the inverse of the ratio 
of the effective area at the photon position and energy to 
the central effective area at the same energy. 

In this paper, we follow the method of background 
subtraction as proposed by Arnaud and collaborators (cf. 
Arnaud et al. 2001; Pratt, Arnaud & Aghanim 2001; 
Arnaud et al. 2002; Majerowicz, Neumann & Reiprich 
2002). Because the cosmic ray background varies from ob- 
servation to observation by ^ 10%, we normalize the back- 
ground level of the blank sky files to that of the source ob- 
servations. The ratios of the total count rates of the source 
observations to those of the blank sky data sets, extracted 
in the whole FOV of each camera in the energy band of 
10-12 keV for the MOS data and 10-15 keV for the pn 
data, are adopted as the background normalization fac- 
tors. They are 1.02, 1.00 and 0.96 for MOSl, M0S2 and 
pn respectively, for this observation. 

The difference in the soft spectral component of the 
cosmic X-ray background between the soiirce region and 
the blank field is investigated as follows. We first accumu- 
late spectra in the annular region of 10-12' for both source 
and blank sky data sets. This region is centered on the 
emission centroid of the group, and is located outside the 
group region. The emission centroid is computed within 
a 2.5' radius. Then we subtract the normalized blank sky 
spectrum from the source spectrum to obtain a residual 
background spectrum, which is subsequently subtracted 
from the spectrum of the group after rescaling its level ac- 
cording to the size of the area where the group spectrum 
is extracted. 

We generate the on-axis ancillary response files 
(ARFs) using SAS task arjgen, and adopt the redistri- 
bution matrix files (RMFs) mll_r7.im_all_2002-ll-07.rmf. 
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Right ascension 

Fig. 1. Vignetting corrected MOS image of RGH 80 in the 
0.5-3.0 keV energy band, adaptively smoothed with a signal 
to noise ratio of 10 and a maximum gaussian smoothing size 
of (7 = 5 pixels. 



in21_r7_im_alL2002-ll-07.rmf and epn_ef20_sdY9.rmf for 
MOSl, M0S2 and pn respectively. 

Following the standard method ^, we have corrected 
the out-of-time events of the pn data by a factor of 0.0232. 

We show the vignetting corrected 0.5-3.0 keV image 
in Figure ^ which is produced by adaptively smoothing 
the combined exposure-corrected MOSl and M0S2 image. 
We find that the X-ray emission of the group is extended 
and almost symmetric, except that it is contaminated by 
several fairly bright sources at 3.5-6.5' from the group 
center. We then use SAS task edetect_chain to perform 
source detections in five energy bands (i.e., 0.2-0.5, 0.5- 
2.0, 2.0-4.5, 4.5-7.5, and 7.5-12.0 keV) and determine the 
source sizes following Katayama et al. (2002), which is 
described below. We divide the whole FOV into 8 annular 
regions centered at the centroid of the group emission, and 
extract an image for each region. The images are binned 
with a bin size of 40 pixels, and smoothed with a maximum 
gaussian smoothing size of cr = 5 pixels. For each annular 
region, the average count rate per pixel and its error (cr) 
are found by recursively excluding the pixels with count 
rate 4(T above the mean value from the image until the 
mean count rate per pixel converges at a constant. Finally, 
for each source detected by edetect_chain except for the 
group, we determine its size, within which the count rate 
is 4tT higher than the mean value. The size of the point 
sources determined by this method is ^ 20". We mask out 
all the sources in our further analysis for both source and 
blank sky data sets. 

see http://wave.xray.mpe.mpg. de/xmm/cookbook/EPIC!PN71 
ootevents.html 



Table 1. Quality of deprojected spectral fits 
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0.58-1.67 
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237 


320 
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4.08-7.67 
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371 





3. Spectral analysis 

While the diffuse X-ray emission of RGH 80 extends 
to a radius of 4.16' at 3a level, the spectra in region 
4.16' < r < 7.67' still show strong features of a ~ 0.95 
keV plasma. We therefore include this region in the spec- 
tral analysis. We extract spectra for a series of annuli cen- 
tered at the centroid of the group emission in such a way 
that for each annulus, with a width larger than 30" (PSF 
consideration), the MOS spectrum has at least 3000 source 
counts in the 0.4-5.0 keV band to give a robust constraint 
on elemental abundances and gas temperatures. Finally, 
we obtain a total of four radially binned spectra for each 
camera (see Table 1 for details) . The contribution from the 
source to the total count rate ranges from ~ 97% for the 
innermost spectrum to ~ 20% for the outermost one. The 
spectra are then grouped so that each bin has at least 20 
counts, thereby allowing statistics to be used. We use 
XSPEC version 11.2.0 (Arnaud 1996) to analyze the spec- 
tra with one or two photoelectrically absorbed VAPEC 
spectral components (Smith et al. 2001) by fixing the hy- 
drogen column density at the Galactic value (1.05 x 10^° 
atom cm^^; Dickey & Lockman 1990). For the spectra 
of the central region, more complicated spectral models 
such as cooling fiow models are also examined. We adopt 
the solar abundances from Anders & Grevesse (1989) with 
Fe/H= 4.68 x 10~^ by number, which is 1.45 times larger 
than the meteoritic value. We have taken this difference 
into account when comparing with other studies. We di- 
vide the elements into six groups, i.e., O and Ne; Mg; Si; 
S and Ar; Fe and Ni; and the others (see also Finoguenov, 
Arnaud & David 2001). In group 1-5, the metal abun- 
dances in the same group are tied to each other and are 
set free in the spectral fittings. In group 6, metal abun- 
dances are fixed at the solar values. We perform a joint fit 
to the MOSl, M0S2 and pn spectra with the same model 
and model parameters, only allowing the normalization of 
each spectrum to vary independently. 

3.1. Deprojection 

With the standard "onion peeling" technique, we calcu- 
late the deprojected spectra by subtracting the emission 
projected from the outer shells, assuming that the spectral 
shape within each shell is the same and that the surface 
brightness profile is described by a /3 model with /3 = 0.45 
and Tc = 0.2', as will be shown in S 14.11 In order to take 
the contribution from beyond the outermost shell into ac- 
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Fig. 2. The deprojected (diamonds) and projected (crosses) temperature and abundance profiles and la errors for the IT 
(black) and 2T (gray) models. The emission measures of the hot and cool components of the 2T model are displayed in Figure|S| 



count, we assume that the spectrum therein has the same 
shape as the one in the outmost shell and has the inten- 
sity proportional to the surface brightness of that region 
(e.g. Matsushita et al. 2002). Note that the deprojection 
analysis will propagate the errors of outer shells to in- 
ner ones. Therefore, the errors of deprojected spectra will 
become larger than those of the corresponding projected 
ones. The reduced chi-squared values of the deprojected 
spectral fitting become smaller than unity in most cases 
accordingly (see Table 1). In what follows We will mainly 
focus on the analysis of the deprojected spectra to obtain 
the three dimensional X-ray properties of the group. 

3.2. Radial temperature & abundance profile 
3.2.1. Single temperature fits 

First we examine the deprojected spectra with a simple 
single-temperature VAPEC model (IT). The best-fit re- 
sults are summarized in Figure[21with black diamonds and 
in Table 1. We find that the temperature increases from 
~ 0.88 keV at the center to - 1.18 keV at 0.6-4', and 
then decreases to ^ 0.95 keV at large radii. Meanwhile, 



all the metal abundances decrease monotonically with ra- 
dius, expect for the Si abundance in the outermost re- 
gion that may have been biased by contamination of the 
background. In particular, the profile of the S abundance 
shows a steeper gradient: From the center to the outer 
regions, Zg/Zpc drops from 2.10 solar to 0.46 solar. This 
relatively higher enrichment of S around the center was 
also found in M 87 with XMM-Newton (Matsushita et al. 
2003) and in the group NGC 1550 with Chandra (Sun 
et al. 2003). Within 7.67', the emission weighted tem- 
perature is 1.01 ± 0.01 keV, and the average iron abun- 
dance is 0.20 ± O.OIZ©, which are in good agreement 
with the previous ASCA results (fcT = 1.02 ± 0.05 keV, 

0.05 ' 
0.06' 



Fe= 0.20tnilZQ; Davis et al. 1999) 



It can be seen from Table 1 that the fitting of the best- 
fit IT model is good for shells 2-4 while becomes much 
worse for the innermost region. As shown in Figure 13 the 
poor fit is due to the large residuals at 0.7-1.2 keV and at 
above 2.5 keV where the model predictions are apparently 
below the observed data. The former is mainly caused by 
the excess Fe L emissions, which implies the existence of 
cooler gas, and the latter indicates that there should be 
a harder spectral component. This suggests the need for 
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a multiphase IGM model (cf. Buote 2000, 2002; Buote et 
al. 2003a). 

We have also made an attempt to fit the projected 
spectra with the IT model and plotted the results in 
Figure 13 with black crosses. In terms of the derived gas 
temperature and metal abundances the difference between 
the projected and deprojected results is only minor, with 
the deprojected analysis giving slightly larger abundances 
(~ 10%) for most cases. 

3.2.2. Two temperature fits 

In order to improve the fits to the observed spectra, we 
employ a multiphase gas model by adding another thermal 
spectral component to the IT model. In this two temper- 
ature model (2T), the two VAPEC components are sub- 
jected to a common absorption that is fixed at the Galactic 
value. The metal abundances of the two components are 
tied to each other, while the gas temperature and normal- 
ization of the second thermal component are left free. 

As shown in Table 1, the 2T model gives significantly 
better fits to the data than the IT model does, especially 
for the inner two shells (e.g. FigureOl). For shell 3, the 2T 
model improves the fitting only slightly. In Figure |21 we 
display the temperature and abundance profiles obtained 
with the 2T model with gray diamonds. The emission mea- 
sures of the hot and cool components of the 2T model are 
displayed in Figure El Within 2', the temperature of the 
cool component remains constant at ~ 0.82 keV, which 
is close to the central temperature obtained with the IT 
model fit. Within measurement uncertainties, the temper- 
ature of the hot component is also constant over the whole 
group region. 

The best-fit of the 2T model gives significantly larger 
abundances than the IT model for O, Mg, Si and espe- 
cially Fe. For example, within the central 1.67', the Fe 
abundance obtained with the 2T model is ~ 0.76Zq, 
which is nearly 2 times higher than that obtained with 
the IT model (~ 0.42Zq). Unlike other metal elements, 
the S abundance does not change notably with the new 
model. Indeed, Buote et al. (2003b) also reported a simi- 
lar trend in group NGC 5044. They attributed the insen- 
sitivity to the fact that the blended helium-like S triplet 
(2.45-2.46 keV) is located farther away from Fc L lines 
(~ 1 keV) than other light elements such as O, Mg, and 
Si, so that the determination of S abundance is less af- 
fected by the change of underlying continuum induced by 
the new model. 

We have also applied the 2T model to the projected 
spectra and shown the results in Figure [3 with gray 
crosses. The projected temperature and abundance pro- 
files are consistent with the deprojected ones within la 
errors in all shells, except that the deprojected tempera- 
ture of the cool component of shell 3 is higher than the 
projected one by ^ 20%. This may arise from the fact that 
the projected contribution of the cool component of the 
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Fig. 4. Detailed single temperature structure. Solid and 
dashed lines stand for the best-fit models to the deprojected 
(solid diamonds) and projected (dotted diamonds) tempera- 
ture profiles, respectively. 

outer shells is removed from the deprojected spectrum of 
shell 3. 

Because in the deprojected analysis uncertainties of 
the abundance determinations are large, and because the 
shapes of the radial abundance profiles obtained with the 
projected and deprojected models are similar, in what fol- 
lows we will use the metal abundances obtained in the 
projected analysis to explore the abundance ratios and 
their constraints on the supernova enrichment scenarios. 
It is found that the IT and 2T models give consistent 
abundance ratios for 0/Fe, Mg/Fe and Si/Fe within Icr 
errors. Moreover, these abundance ratios show no statisti- 
cally significant variations with radius. The derived Si/Fe 
ratios are approximately constant at ~ 1.25 solar for the 
IT model and at ~ 1.05 solar for the 2T model, respec- 
tively. The Mg/Fe ratios are around 1 solar for both mod- 
els. On the other hand, the 0/Fe ratios are subsolar. For 
example, in shell 1 the 0/Fe ratios are ~ 0.58 and ^ 0.55 
solar for the IT and 2T models, respectively. By contrast, 
the S/Fe ratios are different between the two models. As 
it has already been shown in § 13.2.11 the S/Fe ratios ob- 
tained with the IT model decrease significantly with ra- 
dius, while the S/Fe ratios obtained with the projected 
2T model are constant at ^ 0.83 solar within 1.67', and 
decrease slightly in the outer regions. 

3.2.3. Detailed temperature structure of the IGM 

In § 13.2.11 and § 13.2.21 in order to obtain tight constraints 
on spatial abundance distributions, we divide the group 
into relatively few annular regions to avoid poor statistics 
in the spectral analysis, at the cost of losing detailed infor- 
mation on the temperature gradient in each shell. In this 
section, we accumulate spectra for a set of narrower an- 
nuli within 7.67' (462/i^q^ kpc) to perform a detailed study 
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of the temperature structures. After deprojection, we fit 
the spectra with a single temperature VAPEC model, with 
the abundances of O, Mg, Si, S and Fe fixed at the best-fit 
values obtained in S 13.2. f l or their linear interpolation. We 
show the resulting temperature profile in Figure ^ with 
solid diamonds. For comparison, we also show the results 
of the projected analysis with dotted diamonds. 

We find that the temperature profile can be described 
by an empirical lognormal function: 



Table 2. Spectral fits within the central f .67'. 



T{r) = To 



A 
r/ra 



exp ^ 



(Inr — Inro)^ 



(1) 



in which the best-fit values are: Tq = 0.84 ± O.Of keV, 
A = 3.03 ± 0.14 keV, tq = 2.87 ± 0.08', uj = 0.87 ± 0.04 
(;^;2^2.0/2) and To = 0.93±0.01 keV, A = 1.13±0.04 keV, 
ro = 2mtom, ^ = 0.56 ±0.02 (x^=0.8/2) for the depro- 
jected and projected temperature profiles, respectively. 



3.3. Multiphase IGM in the innermost region 



According to our results in § 13.2.11 and § 13.2.21 the 
cool component contributes up to ~ 75% and ~ 30% 
to the 0.4-5.0 keV luminosity emitted by the gas in 
shell 1 and shell 2, respectively. Therefore, these re- 
gions may be important for studying the properties of 
the multiphase IGM as well as investigating the evi- 
dence for cooling flows. To this end, we extract spec- 
tra from the r < 1.67' region, and subtract the con- 
tributions from the outer shells with the method that 
has been described in § 13.11 We then apply several more 
sophisticated spectral models to the data in the 0.4- 
7.0 keV energy band: Model A is a two-temperature 
thermal model, described as WABS(VAPEC-h VAPEC) 
in XSPEC. Next we replace one of the thermal 
components in Model A with an isobaric multi- 
phase component multiplied by an absorber located 
at the source [WABS(VAPEC+ZWABS*VMCFLOW) in 
XSPEC; model B] to examine if the spectra are consis- 
tent with a multiphase cooling flow. Model C is defined as 



Model 


A 


ga 


C 






440.6 


438.0 


450.9 


557.9 
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441 


441 


442 
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Th (keV) 


1 CI +0.07 
J^-OJ^-0.06 


1 03+"°^ 

J^-*^ '5-0. 01 


1 61+°-"^ 


0.98tlt 


Tc (keV) 


0.82lHi 


0.1 (fix) 


tSfi+0 02 




EM^ 


Q (-\fy-\-0.22 

•5-^ ' -0.25 


4 fiq+O-42 
4.D.3_0.43 


r O-1+0.67 

'^•°-'--0.47 


EMi or M"- 


9 qc+O.lS 

Z,.00_QQQ 


on q+3.9 

3 9 


14 6+"® 




O (Zo) 




21+0 07 


26+0 09 


S2+012 


Mg (Zq) 


0.85i°:i« 


11+007 

^^■^^^^-0.07 


0.79l°:°« 


56+0 " 


Si [Zq) 


U.80_o.09 


o.sotro^ 




58+0-" 


S (Zq) 


'^■'■-'-0.20 


0.40t«;i? 


70+0 12 

^- '^-0.21 


0.40l°-J» 


Fe (Zo) 


76+0 08 


0.34lH^ 


70+0 05 


0.48l°:J« 



" The best- fit absorption is O.Sll^ fJ^ x 10^^ atom cm'^ 

The power index and the normalization of the powerlaw com- 
ponent are 1.67!^ojg ^^'^ 2.53!llggg photons/keV/cm^/s, re- 
spectively. 

Emission measure of the hot component, expressed as 
10-'"^ n^nnV/ {An Di{l + zf). 

Emission measure of the cool component of model A. 

Mass deposition rate of the cooling flow in unit of MQyr~^, 
for model B and C. 



WABS(VMCFLOW) in XSPEC. We fix the hydrogen col- 
umn density at the Galactic value. In model A and B, the 
metal abundances of the two spectral components are set 
to be equal. The higher temperature of the VMCFLOW 
component in model B is tied to the temperature of the 
thermal VAPEC component, while the lower one is fixed 
at 0.1 keV. In model C, the two temperatures are left 
free. As pointed out by Molendi & Pizzolato (2001), the 
spectrum of a single-phase gas accumulated in the region 
where the gas temperature gradient is large may also ap- 
pear as multiphase. In this case, the gas is not truly mul- 
tiphase, and therefore model C, which can describe the 
minimum and maximum temperature of the single-phase 
gas within the given region, should give a better fit to the 
data than model A and B. We also examine the spectra 
wfth VAPEC+POWERLAW (model D) to check if the 
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higher temperature component is associated with the cen- 
tral radio galaxies. We show the best-fits of these models 
in Table 2. 

It turns out that, model A and B fit the spectra equally 
well and give a slightly better fit (Ax^ ~ 10) to the 
data than model C. However, the abundances obtained by 
model B are only half of the values given by model A. The 
temperature of the hot component in model A is nearly 
twice as high as that of the cool component. Previous stud- 
ies of cooling flow clusters with ASCA data have already 
found that a two temperature spectral model can well fit 
the spectra of cooling flows (e.g. Ikebe et al. 1999). In par- 
ticular, Ikebe (2001) showed statistically that the ratios 
between the hot and cool temperatures are constant at 
^ 2 for those clusters who demonstrate a very strong cool 
component in the ASCA spectra. This was taken to imply 
that the central two-phase gas reflects the gravitational 
potential structure which has two distinct spatial scales: a 
main cluster component and a second small-scale system. 
In terms of the face values alone, the temperatures of the 
two components we obtained do favor Ikebe's conclusion: 
The hot temperature (~1.5 keV) is close to the virial tem- 
perature of the group, and the cool one (^0.8 keV) is close 
to the kinetic temperature of stars in an elliptical galaxy. 

As has been demonstrated by Bohringer et al. (2002), 
a cooling flow with a broad range of temperature would 
result in a quite broad peak for the blend of Fe L-shcU 
lines. One has to either introduce an intrinsic absorption 
(model B) or leave the lower temperature cut-off to be 
determined by the spectral fitting (model C) in order to 
suppress the excess emission of the low temperature gas. 
In our case, model B fits the spectra better than model C 
with a large intrinsic absorption and mass deposition rate. 
But the large intrinsic absorption column depths may im- 
ply that cooling flow model is an incorrect spectral model 
to the data (Molendi & Pizzolato 2001). Moreover, recent 
results from the RGS instrument onboard XMM-Newton 
have provided little evidence for the X-ray emission from 
the gas with temperature below a certain lower limiting 
(Peterson et al. 2001; Tamura et al. 2001; Kaastra et al. 
2001; Xu et al. 2002; SakeUiou et al. 2002). 

If we allow the hydrogen column density to vary freely, 
we find that the fits are improved, with Ax^ = 14 
and 22 for model A and C, respectively. Model C fits 
the data equally well as model A, with the temperature 
varying continously from 1.49 ± 0.06 keV to 0.56 ± 0.02 
keV. However, the resulting hydrogen column densities of 
model A and C are unreasonably large, which are at least 
6 times larger than the Galactic value. It is unlikely that 
there exits such a large amount of intrinsic absorber in 
this group. 

4. Gas density profile 

In this section, we derive the radial density profile of the 
hot gas in the group, assuming that the gas is single-phase 
with a spatial temperature gradient as found in § 13.2.31 
For the inner region of the group, we will also attempt 
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Fig. 5. Upper panel: cooling function, A, profiles of MOS 
(dashed line) and pn (solid line), scaled to the central values. 
Middle panel: combined MOS and pn projected emission mea- 
sure profile. See text for the detailed derivation of this profile. 
The solid line is the best-fit single /3 model with the convolu- 
tion of the XMM-Newton PSF. Lower panel: electron density 
profiles with la errors obtained from the projected emission 
measure (solid line) and from the normalization parameters of 
the deprojected IT spectral fitting in S I.'-{. 2. .'-{I (crosses'). 



to obtain the gas distribution assuming that the IGM is 
composed of cool- and hot-phase gas. 

4.1. Single-phase gas profile 

The single-phase gas density profile can be easily obtained 
from the observed surface brightness profile. We extract 
the surface brightness profiles in the 0.5-3.0 keV energy 
band for both source and blank sky data sets. The vi- 
gnetting correction and background subtraction are per- 
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Fig. 6. Upper panel: Emission-measure profiles of the hot 
(open circles) and cool (filled circles) components. Black and 
gray lines are the best-fit /3 models. Middle panel: Filling fac- 
tor profile of the cool component and the Icr errors (dotted 
lines). Lower panel: electron density profiles of the hot (gray) 
and cool (black) components and the Icr errors (dotted lines). 

formed using the same method as described in § |2 Since 
the coohng function, A, depends significantly on the tem- 
perature and metal abundance in the temperature range 
found in groups and cool clusters (see Figure EJ, we need 
to take the A profile into account when deriving the gas 
density profile from the observed surface brightness (Pratt 
& Arnaud 2003). We calculate the A profile by using an 
absorbed VAPEC model that is convolved with the in- 
strumental response in XSPEC. The parameters in the 
VAPEC model are fixed at the best-fit projected temper- 
ature and abundance profiles, which have been given in 
§ 13.2.31 and § 13.2.11 From the center to the outmost re- 
gion, the cooling function decreases with radius by 40%. 
Diving the observed surface brightness by the A profile 



gives rise to the projected emission measure, which is pro- 
portional to J n^di. Since the resultant projected emission 
measure profiles of MOSl, M0S2 and pn are well consis- 
tent with each other, we only display the averaged result 
in Figure El 

We fit the projected emission measure profile with the 
standard /3 model (Cavaliere & Fusco-Femiano 1976). If 
no correction of the XMM PSF is made, the projected 
emission measure profile is well modeled by a single (3 
model with fi = 0.419 ± 0.001, = 0.195 ± 0.001', and 
EMo = 1.15 ±0.01 X 10-3 cm-5 {xl = 142.1/82). When 
the convolution of the XMM PSF is performed, we find 
that the fit is equally good, with (3 = 0.415 ± 0.001, 
Tc = 0.105 ± 0.001', and EM^ = 2.68 ± 0.01 x lO'^ cm'^ 
(x^ = 142.9/82). Both values of /3 agree nicely with the 
power-law relation of /3 = (0.44 ± 0.06)r°-20±o o3 found 
by Sanderson et al. (2003). Note that the relatively large 

is caused by the anomalous fluctuations of the surface 
brightness profile, which can not be reduced even if the 
projected emission measure profile is fitted with a two- 
component model. 

Finally, we deproject the /3 model (PSF corrected) by 
using the inverse Abell integral to calculate the gas den- 
sity profile (e.g. Sarazin 1986). The resulting electron den- 
sity profile is plotted in Figure El with the central electron 
density n^Q = 3.27 ± 0.01 x 10~^ cm"^. The rieo is calcu- 
lated from the best-fit parameters of the (3 model using the 
Monte Carlo method: We first generate 10^ random dis- 
tributions of the parameters {EMo, and /3) around the 
best-fit values with the standard deviations found above, 
using a normal probability distribution. We then compute 
TT-oo for each set of the random {EMq, and (3). Finally, we 
obtain the mean of the 10^ rieo and its standard deviation. 
Throughout this paper, we use this method to calculate 
the quantities and their errors, such as the filling factor 
and electron densities in § 14.21 the gas entropy, cooling 
time and masses in § O the SN type fractions and mass- 
to-light ratios in § 16.21 in which the measurement errors 
of all the related parameters are taken into account. 

We have also tried to obtain the gas density profile 
from the normalization parameter, K oc n^V , which is 
obtained in the deprojected spectral fittings with the IT 
model in § 13.2.31 The derived electron density is in good 
agreement with that obtained from the surface brightness 
profiles, except for the central region where the former is 
smaller than the latter by ^ 43%. This again illustrates 
the defect of the IT model. 



4.2. Two-phase gas profile 

In this section, we use the best-fit normalization param- 
eters obtained in the deprojected 2T spectral fittings to 
calculate the gas density distributions for the multiphase 
gas in the inner part of the group. We reproduce four 
deprojected spectra within 2.4', and fit them with a two 
temperature VAPEC model. We fix metal abundances to 
the abundance profiles found in § 13.2.21 and assume that 
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the temperatures of the hot and cool phases are spatiaUy 
constant at the best-fit deprojectcd values of 1.56 keV 
and 0.83 keV, respectively. The emission measure of each 
phase is straightforwardly calculated using the normaliza- 
tion parameters in the VAPEC model and is then modeled 
with a single [3 model (see Figure E)), 



EM{R) ^ EMq 



1 



-3/3 



(2) 



where R is the 3-dimensional radius. 

Similarly to Ikebe et al. (1999), we define the volume- 
filling factor of the cool component /(i?) as EMc{R) = 
n^^^{R)f{R), where rtg^c is the gas number density of the 
cool component. Thus, the emission measure of the hot 



component can be expressed as EMh{R) 



f{R)], where rig^h is the gas number density of the hot 
component. Assuming local pressure balance between the 
two phases, ng^c{R)Tc = ng,h{R)Th, we are able to cal- 
culate the filling factor and subsequent radial gas density 
distributions of the hot and cool components. As shown in 
Figure El f{R) becomes smaller than ~ 0.1 outside ~ 60 
kpc (1'), indicating that the cool gas occupies only a small 
fraction of the volume therein. 

5. Gas entropy, cooling time and mass 
distributions 

We will derive other dynamical quantities of the IGM such 
as gas entropy, cooling time, gas mass and total mass, as 
well as gas mass fraction, using the deprojected tempera- 
ture profile shown in Figure^ and the radial distribution 
of the gas density given in § 14.11 and § 14.21 Following 
the convention in the literature (e.g. Ponman, Cannon 
& Navarror 1999), for the single-phase gas we define the 
gas entropy as S* = k^T /v^^^ , and the cooling time as 
icooi = I ""^khA • turns out that the resulting gas en- 
tropy increases monotonically with radius, and does not 
show an isentropic floor at the central region as far as can 
be resolved with the XMM instruments (Figure [TJ . The 
gas cooling time appears to be shorter than the Hubble 
time within 80 kpc, where the gas has been found to 
have a multiphase nature 13.3(1 . For the two- phase gas, 
we calculate the entropy as 



S = exp 



/nc,c In 5'c + (1 - /)nc,h In 5'h 
/nc,c + (1 - .f)nc,h 



(3) 



where the Sc and are the gas entropies of the cool 
and hot components, respectively. As seen in FigureEl the 
derived entropy of the two-phase gas is larger than that of 
the single-phase gas by a factor of ~ 3 at the center and 
of 1.3-2 for R >20 kpc. 

If the IGM is single-phase, we calculate the gas mass 
by integrating the gas density over the group volume as 
/ nm-prigdV . If the IGM is two-phase, we calculate the gas 
mass of the cool and hot component using J /^TOpng^c/dV 
and J /Ltmpng_h(l — /)dy, respectively, where /i = 0.6 is 
the mean molecular weight, and mp is the photon mass. 




Radius (kpc) 

Fig. 7. Gas entropy (upper) and cooling time (lower) profiles 
with Icr errors for the single-phase gas (black region) and two- 
phase gas (gray region). The solid line represents S{r) oc r^'^ , 
predicted by purely shock heating (Tozzi & Norman 2001). 
The vertical dotted line marks 0.1r2oo, and the horizontal line 
marks the Hubble time. 

The total gravitational mass distribtuion is derived un- 
der the assumption of hydrostatic equilibrium: 



M{< R) 



R^ dP 
'pMdR' 



(4) 



where, P is the gas pressure that is expressed as P = 
n^k^T for a single-phase gas, or P = n^.^k^T = ng i^ksT 
for a two-phase gas. pg is the gas mass density, defined 
as Pg = /iTOpng for a single-phase gas and as pg = 
/iTOp[/ng c -t- (1 — /)ng^h] for a two-phase gas. 

In Figure IHl we display the spatial distributions of the 
gas mass and total gravitating mass, as well as the dis- 
tribution of the gas mass fraction. In 30-80 kpc, the gas 
mass, total mass and gas mass fraction obtained under 
the two-phase assumption are consistent with those given 
by assuming that the gas is single-phase. In the inner re- 
gions, however, the masses derived from the two-phase 
assumption are 3-5 times lower than those obtained with 
the single-phase assumption. This can be attributed to the 
enhanced emission of the cool component in the 2T model. 
Using the single-phase assumption, we find that the gas 



-3/2 



at 



mass fraction rises rapidly with radius to 0.15h^^ 
R - 462/1^0^ kpc. 

We have also examined the effect of the temperature 
gradient on the determination of total mass. We find that 
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Fig. 8. Radial profiles of gas mass, total mass and gas mass 
fraction of the single-phase (black) and two-phase (gray) gas 
with la errors. Dotted lines represent the results obtained by 
assuming the gas is spatially isothermal. 



the isothermal assumption leads to an overestimation of 
the total mass by 17%-30% within 80 kpc, where the gra- 
dient in temperature profile appears to be positive. In 80- 
430 kpc, where the temperature demonstrates a negative 
gradient, the isothermal assumption can result in an un- 
derestimation of the total mass by 48% at most. If the 
observed data are further extrapolated outside 430 kpc, 
the isothermal assumption results in an overprediction of 
the total mass again. At i? ~ SOO/i^q^ kpc, the typical 
value of the virial radius of galaxy group at 1 keV, the 
overestimation is about 17%. 

We have investigated the mass distribution models 
that are suggested by high resolution N-body simula- 
tions, such as the generahzed NFW models (GNFW), 
p{r) = (r/r,)°(i+r/r,)^-° ' ^^'^ ^"^^^ parameters and 



Ts are determined by fitting the observed deprojected tem- 
perature profile with that predicted by assuming that the 
gas (described as (3 model in § I4.1|l is in hydrostatic equi- 
librium with the underlying dark matter-dominated grav- 
itational potential well. Unfortunately, we fail to find an 
acceptable fit with any possible GNFW models by setting 
a = 1 (Navarro, Frenk & White 1997, NFW), or a = 1.5 
(Moore et al. 1999), or a = 2. On one hand, we have only 
6 data points which may be insufficient to give a stringent 
constraint on the models. On the other hand, the tempera- 
ture distribution of this group is complex, for example, the 
existence of two-phase components at the central region, 
so the employment of a single-temperature assumption is 
over-simple. 

6. Discussion 

6.1. Gas temperature profile 

Extrapolating the observed temperature and density pro- 
files to large radii, we can estimate the virial radius 
of the group, r2oo, the radius within which the mean 
mass overdensity is 200 times the critical density of 
the universe. This yields r2oo = 786ft.^Q^ kpc. Such a 
value is, nevertheless, 27% smaller than the one (r2oo = 

1083/i^Q^ kpc) derived from the empirical relation r2oo = 
1.138(T/keV)i/2(i _|. 2)-3/2/j-i j^^^^ j^Q^ ^j^g 

merical simulations (Navarro, Frenk & White 1995). 
Regardless of the fact that this empirical relation may 
significantly overpredict r2oo in small halos (Sanderson et 
al. 2003), we will adopt r2oo = 1083/i^o^ kpc below, un- 
less stated otherwise, for the convenience of comparing 
with other studies. In this case, the observed diffuse X- 
ray emission of the group extends to ~ 45% of the virial 
radius. 

The temperature profile of RGH 80 follows more or 
less a universal form as has already been found in some 
groups of galaxies observed with ROSAT (Mulchaey 2000, 
and references therein) and ASCA (Finoguenov et al. 
2002a): Gas temperature is at its minimum (0.83 keV) at 
the center, rises to the maximum value at a small radius 
(0.11r2oo), and then drops gradually with radius. Similar 
temperature profiles have been found recently with XMM- 
Newton and Chandra in galaxy groups NGC 1550, NGC 
2563, NGC 5044, and MKW 4 (Sun et al. 2003; Mushotzky 
et al. 2003; Buote et al. 2003a; O'Sullivan et al. 2003). 
Yet, we note that there are still some groups in which a 
fiat temperature profile at large radii is reported, such as 
NGC 1399, NGC 2300, NGC 4325 and AWM 4 (Route 
2002; Mushotzky et al. 2003; O'Sullivan & Vrtilek 2003; 
see Sun et al. 2003 for a comprehensive summary). 

6.2. Gas metallicity and supernovae enrichment 

With the abundance ratios derived in § 13.2.11 and § 13.2.21 

we are able to probe the origin of the heavy elements in 
RGH 80 and to distinguish between the different contribu- 
tions of type la and type II supernovae. Following Gibson, 
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Fig. 9. The radial profiles of SN type fractions derived from the observed abundance ratio pattern. Only the results of those 
models in which consistent SN type fractions can be achieved with the observed Si/Fe, S/Fe, O/Fe and Mg/Fe ratios are shown. 



Loewenstein & Mushotzky (1997; sec also Ishimaru & 
Arimoto 1997), we take MFe,sNia/-^^Fc, total to quantify the 
roles played by SN la in contributing to the metal abun- 
dances. The relative frequency of SN la to the total SNe, 
C, can be determined by the ratio of the total mass of the 
a elements to the total iron mass, together with the SN 
models. 

For SN la, we adopt seven updated models listed in 
table 3 of Iwamoto et al. (1999), namely, the W70 model 
in which the initial mctallicity is assumed to be zero, the 
W7 model and five delayed detonation models, WDDl-3 
and CDDl-2. We also utilize the SN II model listed in the 
same table (the 199 model hereafter), which is similar to 
the T95 model in Gibson et al. (1997; see also Tsujimoto et 
al. 1995). In addition, we also choose the four representa- 
tive "W95" models for SN II (W95;A;1O-4Z0, W95;A;Z0, 
W95;B;10'4Zq, W95;B;Zq; Gibson et al. 1997; Woosley 
& Weaver 1995). For all the SN II models, the elemental 
yields are integrated over 10 Mq to 50 Mq by using the 
Salpeter IMF. 

We first utilize the metal abundances obtained from 

the spectral fittings with the projected 2T model, in which 
the abundances are better constrained than those with 
the IT model at the central region. Based on the best- 
measured Si/Fe ratio, we find that 90% 100% of the iron 
mass was produced by SN la according to the SN la model 



WDDl/CDDl together with any of the SN II models. The 
fraction decreases to 70%-90% if the other combined SN 
models are used. The SN la fractions estimated with the 
observed O/Fe ratio agree well with those obtained with 
the Si/Fe ratio. The fractions range from 65% to 90% for 
any combinations of the SN models. 

The S/Fe ratio, on the other hand, gives a somewhat 
loose constraint on the SN la contribution, which ranges 
from 65% to 100% for any possible combinations except 
for those combined with WDDl and CDDl. These two 
SNIa models (WDDl and CDDl) seem to overproduce S 
relatively to Fe, which results in a unreasonably high SN 
la fraction (> 1). With the observed Mg/Fe ratio and the 
199 model, we obtain a consistent SN la fraction with that 
derived from the O/Fe, S/Fe and Si/Fe ratios. However, 
if the four W95 models are in use, we tend to underesti- 
mate the Mg yields, which leads to a relatively low SN la 
fraction. 

Finally, we find that in those combined models where 
the SN la model is not WDDl or CDDl, consistent SN 
type fractions can be achieved with the observed abun- 
dance ratios of Si/Fe, S/Fe and O/Fe. If the model combi- 
nation is 199-hW70, I99-hW7, I99-hWDD2, 199-hWDD3, 
or I99-I-CDD2, consistent SN la fractions can be ob- 
tained with the abimdance ratios of Si/Fe, S/Fe, O/Fe 
and Mg/Fe. If the model combination is W95;B;ZQ-hW70, 
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W95;B;Z0+W7, or W95;B;Z0+WDD3, a marginally con- 
sistent SN la fraction can be obtained in i? < 20 kpc with 
the above four abundance ratios (see Figure 13 ■ 

In terms of the radial behaviors of the SN la fractions 
derived from the abundance ratios shown in Figure El 
it appears that the SN la fraction remains almost con- 
stant with radius at ^ 85% for the SNII model 199 or at 
- 72% for the SNII model W95;B;Z0. It is likely that the 
decline in the SN la fraction with increasing radius de- 
rived from the Mg/Fe ratio is an artifact because Mg lines 
are blended with the outskirt of the Fe L-shell complex, 
and therefore are measured less accurately than O and 
Si. The resulting SN type fraction agrees with those ob- 
tained from previous ASCA observation of galaxy groups 
(Finoguenov & Ponman 1999) and from recent XMM- 
Newton and Chandra observations of NGC 5044 (Buote et 
al. 2003b), NGC 1399 (Buote 2002) and NGC 1550 (Sun 
et al. 2003). 

We then employ the abundance ratios obtained in the 
projected IT model fittings to examine the SN models and 
contributions to the iron mass of different types of SNe. 
We find that only in models I99-f-WDDl and I99+CDD1 
can a marginally consistent SN la fraction of ^ 80% be de- 
rived from the observed abundance ratios of Si/Fe, S/Fe, 
0/Fe and Mg/Fe. 

We can calculated the ratios of the metal mass in the 
IGM to the total blue luminosity of the galaxies in the 
group for the measured elements within the observed re- 
gion of 462/i5o^ kpc using the elemental abundances ob- 
tained from the projected IT model fittings. The metal 
mass of element i is calculated by 

M,{<R)= / ATrm,ZQ,,Z,{r)nii{rydr, (5) 
Jo 

where m.i is the atomic mass of element i. For the to- 
tal blue luminosity, we take the magnitudes rrizwi in the 
Zwicky B(0) system from Ramella et al. 1995, and convert 
them to magnitude B with the relation mz„i — B + 0.35 
given by Gaztanaga & Dalton (2000). The metal mass-to- 
light ratios within the virial radius of kpc are also 
derived by extrapolating the gas distribution to r2oo and 
assuming that the abundances beyond the observed region 
are the same as those in the outmost shell. We summarize 
the results in Table 3. 

It turns out that the derived Fe M/L is ~ 2 times lower 
than the corresponding typical value of clusters (Arnaud 
et al. 1992; Loewenstein & Mushotzky 1996; Renzini 1997; 
Finoguenov et al. 2001). Moreover, Si M/L between clus- 
ters and this group differs by a factor of ~ 10. This is 
consistent with the findings of Finoguenov et al. (2001). 
The lower values of metal mass-to- light ratios of the group 
may imply that the group has lost some of the enriched 
gas that was produced by the galaxies in the group. 

6.3. Entropy 

The observed gas entropy of RGH 80 at 0.1r2oo is ~ 
90h^Q^^ keV cm for a single-phase gas model, which 



seems to be a typical value for groups of '--^ 1 kcV (Lloyd- 
Davis et al. 2000; Ponman, Sanderson & Finoguenov 
2003). The entropy reaches ^ ISbh^^^^ keV cm if gas 
is assumed to be two-phase. Additionally, gas entropy in- 
creases with radius as r'^-^^ in the region of r > O.Olrvir 
and as r^" in the region of 30 kpc< r < 100 kpc for 
the single-phase and two-phase assumptions, respectively. 
As a result, the derived radial entropy profiles from the 
single-temperature assumption are flatter than that ex- 
pected from purely shock heating (S'(r) oc r^'^; Tozzi & 
Norman 2001) throughout the whole group. This implies 
that the IGM may have suffered from non-gravitational 
effects such as radiative cooling, heating by supernovae 
and/or AGNs, even at large radii. The preheating sce- 
nario, in which gas is pre-heated and subsequently col- 
lapses adiabaticly, predicts an isentropic core of 50-100 
keV cm^ within O.lrvir for groups of 2 x 10"'^'^ Mq in order 
to reconcile the theoretically expected Lx-T relation with 
observation (Tozzi & Norman 2001). However, the entropy 
of RGH 80 we obtained starts to increase rapidly from an 
inner radius of as small as 0.01-0.02rvir, which disagrees 
with the prediction of preheating model. 

The expected entropy profile based on radiative cool- 
ing model or internal heating model, or cooling plus in- 
ternal heating model (Finoguenov et al. 2002b; Borgani 
et al. 2002; Brighenti & Mathews 2001) resembles the 
observed data of this group in shape. In particular, it 
has already been found that the cooling may marginally 
reproduce the observed scaling relations of S{0.1ryir)-T 
(Voit & Bryan 2001) and L^-T (Muanwong et al. 2002; 
Wu & Xue 2002a), and may even be responsible for the 
scale-dependence of the IGM mass function (Wu & Xue 
2002b). However, the cooling process suffers from the so- 
called cooling crisis (Balogh et al. 2001), and is also ineffi- 
cient in the explanation of the observed X-ray properties 
of groups and clusters (Bower et al. 2001). Inclusion of a 
suitable energy feedback from galaxy formation has been 
suggested to prevent the IGM from overcooling and sup- 
ply the IGM with additional energy (Borgani et al. 2002). 
Indeed, Xue & Wu (2003) pointed out that a combination 
of cooling and heating can explain simultaneously the ob- 
served global X-ray properties of groups and clusters (e.g. 
gas entropy distribution and Lx~T relation), and the ob- 
servational limits on the contribution of the diffuse IGM in 
virialized halos to the X-ray background within the frame- 
work of standard CDM structure formation with an am- 
plitude of matter power spectrum as ~ 0.7. Therefore, 
our observation seems to favour the heating plus cooling 
model. 

6.4. comparison with scaling relations 

We now compare the total mass and gas fraction of 
RGH 80 with the expectations from the scaling relations 
of /gas-r and M-T derived by Sanderson et al. (2003), 
adopting a virial radius of 786/i^q'^ kpc. The gas fractions at 
0.3r2oo and r2oo are O.OGh^^ and 0.25/i5Q , respectively. 
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Table 3. The IGM metal mass to light ratios 



Radius 
(kpc) 


Mgas 


Lb 
(1O"L0) 


(W-^Mq/Lq) 


Mo /Lb 
(W-^Mq/Lq) 


Mms/Lb 
(IQ-^Mq/Lq) 


Msi/Lb 
(IO-^Mq/Lq) 


Ms /Lb 
{10-^Mq/Lq) 


462 


3.10 ±0.01 


2.89 


3.75 ± 0.24 


1.19 ±0.75 


1.22 ± 0.45 


1.94 ± 0.50 


9.20 ± 5.33 


786 


7.88 ± 0.03 


4.39 


4.67 ±0.38 


1.91 ±1.56 


1.10 ±0.66 


3.00 ± 1.03 


13.95 ± 11.01 



which are comparable to the upper limits of group gas frac- 
tions of in the Sanderson et al. (2003) sample. This indi- 
cates that RGH 80 is a relatively gas rich system. The total 
masses of the group (1.53 x lO^^h'^^MQ within 0.3r2oo and 
3.14 X lO^'^h^Q Mq within r2oo) agree nicely with what arc 
expected from the observed M{0.3r2oo)-T and M{r2oo)-T 
relations of groups and clusters (Sanderson et al. 2003), 
respectively. 

The velocity dispersion and bolometric X-ray luminos- 
ity of RGH 80 are 450 km s'^ (Ramella et al. 2002) and 
2.06 X lO^^ft-gQ^ erg s~^ within 4:62h'^^ kpc, respectively. 
These values are comparable to the expectations from the 
observed Lx-<7 relation for groups and clusters (Wu, Xue 
& Fang 1999), though the group appears to be relatively 
luminous in the Lx-T plane (Xue & Wu 2000). 

7. Conclusions 

We summarize below the main conclusions from our anal- 
yses of the XMM-Newton observations of the galaxy group 
of RGH 80. 

— The X-ray emission of the group is detected out to 
~ 7.67', or 462/i^g^ kpc, corresponding to 0.45r2oo- The 
group seems to be relatively gas rich and luminous with 
respect to its temperature (1.01 ± 0.01 keV). 

— Spectral analysis shows that the temperature profile, 
which increases from the center and then declines with 
radius after reaching a plateau around 0.11r2oo, fol- 
lows a universal profile (Mulchaey 2000) and that the 
abundance profile of each measured element decreases 
monotonically with radius. 

— In the central region, the X-ray emission of the gas is 
better modeled by a two-temperature spectral model, 
with temperatures of 0.82 and 1.51 keV and the 
Galactic absorption, than by a single temperature or 
cooling flow models. Beyond ^ 60 kpc (1'), the vol- 
ume filling factor of the cool component, f{R), be- 
comes smaller than ~ 0.1, indicating that the cool gas 
occupies only a small fraction of the volume therein. 

— The gas entropy distribution derived from single- 
temperature assumption deviates from the prediction 
of the over the whole observed region, The isentropic 
core at the center expected from preheating model does 
not show up. Nevertheless, our derived entropy pro- 
file resembles what is predicted by radiative cooling 
model or internal heating model or cooling plus heat- 
ing model. 

— Both the abundance ratios of Fe/0 and Fc/Si are very 
high showing a large SN la dominance, which is higher 



than that for M87 (Matsushita et al. 2004) and compa- 
rable to or slightly higher than that for Ccntaurus clus- 
ter (Matsushita et al. 2004). With the abundance ratio 
pattern of the hot gas, we estimate that ^ 85% (199) or 
^ 72% (W95;B;Z0) of the iron mass is contributed by 
SN la. This SN type fraction remains almost constant 
against radius. We find that S is significantly overpro- 
duced by the two delayed detonation models WDDl 
and CDDl, while Mg is slightly underproduced by the 
four W95 models considered here. 
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